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1. Introduction 

The renormahzation group (RG) theory is suitable to the study of many modern physical 
problems. Generically, every situation where the scale of typical physical interest belongs 
to a (wide) range of correlated or coupled scales may be (must be?) treated by RG 
techniques. Critical phenomena, which are characterized by one (or several) diverging 
correlation length(s), provide "the" didactic example [T]. Quantum field theory, with 
its strongly correlated quantum fluctuations, is not less famous since it has given rise 
to the early stages of the RG theory [2j. 
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Thanks to a fortunate success (essentially due to an impressive diagrammatic 
calculation in estimating the critical behavior of some systems jHE], the perturbative 
framework has pushed into the background the undoubtedly nonperturbative character 
13 EI of the RG theory. As a consequence, there has been relatively little interest in 
the development of nonperturbative RG techniques In particular the formulation 
of the RG theory via an infinitesimal change of the scale of reference (running scale), 
designated by the generic expression "exact renormalization group equation" pP although 
known since 1971 ^0], has actually been actively considered only since the beginning of 
the nineties jH]. Because the variety of systems to which the exact RG formulation could 
be applied is large [HI [H] (see also |T2|) and also because the perturbative framework 
is generally not well adapted to such studies 5 it is worthwhile making every 

endeavour to master, if possible, the exact RG framework. 

The exact RG equation is an integro-differential equation the study of which calls for 
approximations and/or truncations. Among the possible approximations, those based 
on expansions in powers of a small parameter such as e = ci^ — or 1/A^ (where the 
upper-dimension (i^ = 4 for the A^-vector model) are perturbative in essence. They, 
however, present the advantage of allowing analytic calculations but are attached to the 
smallness of quantities that are actually not small in the cases of physical interest. In 
some cases the perturbative framework may fail ^1 |H] ■ 

The derivative expansion [TF, of present interest here, is an expansion in powers 
of the derivative of the field. It is not associated to a small parameter though it is 
expected to be rather adapted to the study of phenomena at small momenta J (large 
distances) like critical phenomena for instance. The interest of the derivative expansion 
is that the physical parameters (like d and A^) may take on arbitrary values. Hence, 
in the range of validity of the expansion (thus presumed to be in the large distance 
regime), the approach is actually nonperturbative. The drawback is the necessary 
recourse to numerical techniques [for studying coupled nonlinear ordinary differential 
equations (ODE)] that are not always well controlled. Consequently very few orders of 
the derivative expansion have really been explicitly considered. 

Many studies have been effectuated in the local potential approximation (LPA), 
i.e. at the leading order [O (d^)] of the derivative expansion iniinillE]. Also, estimates 
of the critical exponents for Ising-like models (A^ = 1) have been obtained several times 
from full studies of the next-to-leading order § [i.e. O (d^)] US IHl 122 123 121 CS] , and 
even from a full study of the third order [O [d^)] |2S]- On the contrary, only two full 

J A recent interesting attempt to adapt the derivative expansion to phenomena that include effects at 
larger momenta is made in reference |17|. 

§ In incomplete studies some contributions to a given order of the derivative expansion are neglected. 
For example in reference jl8^ , despite an estimation of the critical exponent 77 (which is exactly equal to 
zero at order 9*^), the order has not been completely considered because the evolution equation of the 
wave function renormalization has been neglected. In reference |19) . O (N) systems are incompletely 
studied up to O (9^) because one differential equation has been discarded. Also, the study of reference 
[20] at order 9*, though interesting, is incomplete since only three differential equations, among five 
constituting actually the order 9^, have been treated. 
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studies have been effectuated up to O (9^) for the O (A^) vector model (2nil2Z|- Yet, this 
model provides us with the opportunity of varying and of comparing the results with 
the best estimates of critical exponents obtained from six fl^ EH] or seven jHOI orders 
of the perturbative framework in a wide range of values of A^. Interesting informations 
on the possible convergence of the derivative expansion are then reachable. 

In the present study, I consider the O (N) Wilson-Polchinski exact RG equation 
expanded up to order O [d"^) in the derivative expansion and I calculate the critical 
exponents. Beyond the specific calculation (which was lacking) the real aim is to try 
to clarify (and also to evaluate) the present status of the derivative expansion. Still, I 
must explain why I consider the Wilson-Polchinski equation. 

Indeed, there are several different approaches and treatments of the exact RG 
equation so that it is not easy to really estimate the validity of the choice of the equation 
or/and of the calculations available in the litterature. Let me try to briefly summarize 
the situation and to justify my choice. 

There are two families of exact RG equations (for a review see j^). The first 
family expresses the RG flow (Ad^A/dA) of the microscopic action [0] (hamiltonian) 
associated to a running momentum scale A which, in the circumstances, is a running 
ultra-violet cutoff. The second family expresses the RG flow of Fa [M] , the Legendre 
transform of ^a [(f)], in that case the running scale A effectively appears as an infra-red 
momentum cutoff. 

There is no fundamental difference between the two families since the object of 
the RG is the same in the two cases: accounting for all the correlated scales A ranging 
from to oo (at criticality); only the physical meaning of the field variable at hand 
has been changed: is related to a microscopic description (like a spin of the Ising 
model) while M is thought of as a macroscopic variable (like the magnetization). If one 
wants to calculate an equation of state or some correlation functions or some universal 
critical amplitude ratios, the second family of equations is better adapted. But if one 
only wants to estimate critical exponents (for example to illustrate the convergence 
of the derivative expansion), then considering the first family is surely more efficient. 
Actually, the set of ODE generated in the derivative expansion is much simpler when 
considered with the first family than with the second ||. The first family is indeed 
better adapted to the calculation of critical exponents for the same reason as in the field 
theoretical approach to critical phenomena j3], the critical exponents are defined from 
renormalization functions that are introduced within the microscopic action Sa [(p] ■ The 
Wilson and Polchinski [22] exact RG equations belong to the first family, they only 
differ by the way the smooth cutoff function has been defined (a specific function for 
Wilson, an arbitrary one for Polchinski). Of course, the two equations are physically 
equivalent, but due to a misunderstanding in the introduction of the critical exponent 

II In the study of reference jSH) of the Ising case up to O (9'') , the consideration of an exact RG equation 
of the second family yields a set of ODE the writing of which requires 20 pages [33 while, even at order 
d^, the first family yields a set of equations that holds on a half of page [321 ■ The complexity of dealing 
with the second family is also well illustrated in the appendix of reference |26| . 
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T] in the exact RG equation as formulated by Polchinski, it is only recently that the 
equivalence has been clearly established IT^. 

If the coupled set of ODE generated by the derivative expansion is different in 
the two families of exact RG equations, the treatments of the differential equations 
encountered in the litterature differ also. For convenience, let me classify the studies in 
two groups according to whether the authors have adopted the conventional approach 
(defined below) or not. 

The conventional approach is characterized as follows: 

(i) the set of ODE is numerically studied as such (e.g., without considering any artefact 
such as an expansion in powers of the field). 

(ii) the critical exponent rj is introduced in a conventional way as defined in references 

(iii) the critical exponents are estimated from a set of eigenvalue equations linearized 
about a fixed point solution of the flow equation. 

(iv) the reparametrization invariance is explicitly accounted for. 

To limit myself to the studies mentioned above that consider O (N) systems via 
equations of the second family developped up to 0(9^) [SEIEZI: the study of reference 
|26j follows the conventional approach while that of reference j^Z] does not. In fact, in 
this latter work, except the first point above, none of the other points is satisfied. This is 
particularly important relatively to point (iv) because the reparametrization invariance 
induces a line of equivalent fixed points along with rj is constant 36 . In the case where 
the invariance is broken (it is generally the case within the derivative expansion, except 
in reference ^B]) then the fixed points along the line are no longer equivalents and 
the effective t] (when introduced conventionally) varies with the global normalisation of 
the field 0. Nevertheless, even if the invariance is broken, one expects that a vestige 
of this invariance can still be observed |36t |HZl Ej via an extremum of r) on varying 
the global normalization of the field 0. The absence of explicit consideration of the 
reparametrization invariance in reference j2Z| is intringuing inasmuch as the estimation 
of the critical exponents is excellent (see section 

By considering the Wilson-Polchinski RG equation for 0{N) systems, my first aim 
is to illustrate the conventional treatment as described above. Additionnal aims follow: 



Morris and Turner ^| have imposed the reparametrization invariance by choosing 



a specific cutoff function and the resulting estimates of the critical exponents are not 
very good, especially for u. Does the Wilson-Polchinski RG equation also produce 
such bad results at order O (9^)?. 

(ii) Is the Wilson-Polchinski RG equation able to produce estimates of critical 
exponents comparable to those obtained by Gersdorff and Wetterich |2Zj? 

(iii) The O (A^) exact RG equation expanded up to order involves three coupled ODE. 
Compared to the Ising case = 1, it is somewhat an intermediate state between 
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the order d"^ which involves two coupled equations and the order d'^ which involves 
five coupled equations [23 1201 • 
(iv) Does one observe some signs of convergence of the derivative expansion already at 
order 

Relative to this latter point, it is worth mentioning here the work of Litim 
whose aim, though very interesting, differs from that of the present work. Litim 
focuses its attention almost exclusively on the second family of the exact RG equations 
and especially on the arbitrariness introduced by the regularization process (cut-off 
procedure) but he does not account for the reparametrization invariance. From general 
arguments (independent from the derivative expansion), he provides us with a criterium 
for choosing the regulator which should optimize the convergence of the derivative 
expansion already at a very low order. The considerations are surely useful especially 
for studying sophisticated systems, such as gauge field theory for example, for which 
already the leading order of the derivative expansion is difficult to implement. However, 
Litim has not actually studied the convergence properties of the derivative expansion 
in itself but, in fact, has implicitly assumed that it converges (at least, the fact that 
the expansion could yield only asymptotic series has been excluded). Furthermore, 
Litim's criterium of choice does not apply to the Wilson-Polchinski RG equation. In 
particular, at leading order of the derivative expansion, his optimization provides the 
same estimates of critical exponents as those obtained with the Wilson-Polchinski RG 
equation (see section 13. 2|) which, at this order, does not display any dependence on the 
regularization process. 



2. Derivative expansion up to 0(9^ 



2.1. Flow equations 

According to reference JH], after subtraction of the high temperature fixed point 



satisfied by an O (A^)-symmetric action S [(p"'] (with a 



from the action, the Wilson-Polchinski exact RG equation 



N) reads as follows: 



S 




2^' ^ \ 



wP-q^P']^^ 



6^S 



6S SS 



(1) 



in which S stands for dS/dt = — AdS/dA (hence exp (— t) = A/Aq in which Aq is 
some initial momentum scale of reference [13j), q is a dimensionless (i- vector {d is the 



dimension of the euclidean space and q = 



d}) , q' 



y-d 2 



P(g2 



IS 



a dimensionless cutoff function that decreases rapidly when g — > oo with -P(O) = 1 , 
ip {q^) is an arbitrary function (except the normalization (0) = 1) introduced to test 
the reparametrization invariance, vj = 1 — r]/2 and = d/2 + w . A prime denotes a 
derivative with respect to q^: ijj' = di/j/dq"^, P' = dP/dq^ and q-<9g/ (g) = qidf/dqi. 
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The expansion up to order d"^ consists in projecting equation 
the form: 

vip,t)+z{p,t)id$y+Y{p,t){$d$f 



onto actions of 



with 



P 



1 ^ 

^ OL = \ 



21 



Then the flow equations for V ^ Z and Y read as follows: 

V = Jo {NV'^2pV") ^dV -{d^ 2zu) pV - 2wpV''^ 

+ 2/1 {NZ + 2pY) 
Z = lo {NZ'+2pZ" + 2Y)-2{w+ 1) Z- {d + 2w) pZ' - 2^'qV' 

- Aw {V'Z + pV'Z') - [{w - 1) Po + 2t^^o] 

Y = Io {NY'+2pY"+AY') -{d + 2 + 4w)Y- {d + 2w) pY' - 2^'qV" 

- Atu {V"Z + 2V'Y + pV'Y' + 2pV"Y) 



(2) 



- 2 [(ro - 1) + 2wil)'^ (V"V' + pV"^) 

in which a prime acting on V, Z ot Y denotes this time a derivative with respect to p, 
while = ^' (0), pQ = P' (0) and: 



zuP-q^P')ip 



q'{wP-q^P')ij 



It is convenient to perform the following changes: 
p = NIop, V = NIoV, Z = Z, 



Y 



Y 



(3) 
(4) 



V2 



V3 = F 



and then to consider the new set of functions: 
dV 
dp 

Using these new notations and restoring the writing p 
(j21) becomes: 

Vi = 



p, the set of equations 



1 + — j v;+— pv'/ - evi -{d + e) py[ - e (v^ + 2pviv; 



N 



+ Pi 



V2 + ^ (V3 + PV3) 



2 2 

V2 = V2+— pv2 +— V3 - (e + 2) V2- {d + e) pv2 + 2mvi 



V3 



- 2e (viV2 + PV1V2) + P2V1 

(1 + ^) ^'s+^P^ -{d + 2 + 2e)v3- {d + e) pv^, + 2uy[ 

- 2e [v[ (v2 + 2pv3) + vi (2v3 + pvg)] + 2P2 (vi + pv[) v[ 



(5) 

(6) 
(7) 



in which: 



u 
Pi 



^'0 



0' 



2w 



--l)P'(0) + en 
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2.2. Fixed point equations 



The fixed point equations correspond to the three simultaneous conditions Vj 
i = 1, 2, 3 which yield three coupled nonlinear ODE of second order each: 



for 



N 



evi(l + vi; 



N 



V3 



N 



+ y (c? + e + 2evi)v;-Piv^ 



N 



(e + 2 + 2evi) V2 - V2-— V3 - 2mvi - P2V1 



N 

+ ^[{d + e) +2evi] 



N 

Yp 



(d + 2 + 2e + 4evi) V3 - ( 1 + — 1 v 



(8) 



(9) 



3 + 2 (ev2 - u- P2V1) v'l 



+Tr [(rf + e + 2evi^ 



2(2ev3-P2v'i)v;] 



(10) 



The differential system is of order six, thus the general solution depends on six 
arbitrary constants. Three of these constants are fixed so as to avoid the singularity at 
the origin p = displayed by the equations, hence the three following conditions: 



(0) = [e + 2 + 2evi (0)] V2 (0) --V3 (0) - vi (0) [2u + P2V1 



v'l (0) 
v:,(0) 



N 



N + 2 
N 



evi (0) [1 + vi (0)]-Pi 



(0) + |^V3 (0)] } (12) 



{[rf + 2 + 2e + 4evi(0)]v3(0) 



N + A 

+2[ev2 (0)-^i-P2Vi (0)]v;(0)} 



(13) 



If ?7 is a priori fixed, then the general solution of the set of equations (jHtiTUI) depends 
on the three remaining arbitrary constants, e.g. the values Vj (0) for i = 1,2, 3. In general 
the corresponding solutions are singular at some varying p* (moving singularity), with: 

vi (p) oc (p* - p)"^ , V2 (p) oc (p* - p)"^ , V3 (p) oc (p* - py^ 

But, the equations (IHtiTUj) admit another kind of solution that goes to infinity 
(p 00) without encountering any singularity and which behaves asymptotically for 
large p as follows: 

viasy (p) = Gip'^ + (1 + 29^) + ■■■ 



2 (rf - e) u + {d+e)P^ 29, 

Yd ^ 



V2asy (p) = uGip'^ + Cj (1 

+ G2p'' + ■■■ 

V3.. (P) = u9,G,p^^-^ + 2G19, nd-e)u + dP'y ,,^_, 



2 (d + e) 



+ G,p 



03 



(14) 



(15) 



(16) 
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with: 



9l = -— , 02 = , 6-i = 62-I 

a + e a + e 



The values of the three constants {Gfji = 1,2,3}, correspond to some adjustment 
of the set {vj (0) ; i = 1, 2, 3} and vice versa. This nonsingular solution is the fixed point 
solution which we are interested in. When 77 is a priori fixed, the six arbitrary constants 
of integration are then determined, the differential system is balanced. 

If rj is considered as an unknown parameter to be determined, then one of the three 
preceding quantities {vj (0)} or {Gi} must be promoted to the rank of a fixed parameter 
chosen a priori. In general one chooses 

V2 (0) = Zo (17) 

which corresponds to having fixed to Zq the value of the kinetic term in 5* [0"] and 
is customarily associated with the arbitrary global normalization of the field 0. One 
thus obtains a function rj (Zq) which should be a constant if the reparametrization 
invariance of the exact RG equation was preserved by the derivative expansion presently 
considered {rj should be a constant along a line of equivalent fixed points generated by 
the variation of Zo). Since it is not the case, one actually obtains a nontrivial function 
f] (Zq). Fortunately a vestige of the reparametrisation invariance is preserved and 77 (Zq) 
displays an extremum in Zq. This provides us with an optimal value (?7°p*) of rj (and 
similarly for Zq) via the condition: 

(^«' 

dZo 

Instead of using this condition to determine 77°?*, I use the fact that the line of 
equivalent fixed points is associated with a redundant operator with a zero eigenvalue^ 
|36| EZj. Hence, one may determine 77°^* by imposing that the fixed point of interest 
be associated to a zero eigenvalue. This leads us to the consideration of the system of 
eigenvalue equations. 



2.3. Eigenvalue equations 

The eigenvalue equations are obtained by linearization of the flow equations (jSHZj) about 
a fixed point solution {v*, i = 1,2, 3}: 

Vi = Y* + ee^'gi 

% The fact that the eigenvahie takes on a definite value although it is associated with a redundant 
operator is not in conflict with the work of Wegner (see reference |39j and references therein) which 
indicates that the eigenvalue of a redundant operator generally varies with the renormalization process. 
As shown in reference ^36_, the linear character of the renormalization of the field in the process of 
generating the exact RG equation implies a definite eigenvalue (the reparametrization invariance is 
a direct consequence of this linearity). On the contrary, in the case of a nonlinear renormalization 
scheme, the eigenvalue in question no longer is constant but depends on the renormalization procedure 
[HHj in accordance with P^ . 
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Keeping the linear contribution in e, the following set of coupled ODE comes: 



N 



(A + e + 2evt)gi- 1 + 



ll - A ( g2 + ^g3 



N 



+ -[{d + e + 2eY\) g[ + 2egivn - Pig^ 



N 
2p 



(A + e + 2 + 2eYl) g2 - ga-^ga + 2 (ev^ -u- P^vl) gi 



+_[(rf + e + 2evng2 + 2egivr] 



(19) 



(20) 



2p 



(A + rf + 2 + 2e + 4ev*) 



2ug[ 



+2e (g'.v; + 2giv*) - 2P2 (givf + v^gi) + 2evf ga] 

+ - [id + e + 2evl) g', + 2e i2g[v; + giv^ + 2vrg3) - AP2g[vr] (21) 

Similar considerations as those relative to the determination of the six integration 
constants associated to the fixed point equations (jSHTUj) stand. For a given fixed point 
solution (v*), there remain six constants to be determined. Three of them are fixed so 
as to avoid the singularity at the origin p = displayed by the equations: 



ll (0) 

?3(0) 



[A + e + 2 + 2evt (0)] g2 (0) --ga (0) 
+ 2[ev; (0)-n-P2vI (0)] gi (0) 
^ { [A + e + 2evl (0)] gi (0) - Pi [g^ (0) + |^g3 (0) 



A^ + 2 
A^ 



{[X + d + 2 + 2e + AcyI (0)] ga (0) - 2ug[ (0) 



A^ + 4 

+ 2e[g'i(0)v; (0) + 2gi (0)v* 
-2P2 [gi (0) vf (0) + vt (0) g[ (0)] + 2evr (0) g2 (0)} 

One is interested in the solution that is regular when p — > oo. 
For A a priori fixed, the three values {gi (0) , i = 1, 2, 3} at the origin p = must be 
adjusted so that the solution reaches the following regular asymptotic behavior: 



with: 



?lasy = S,p^^ +2(1 + ^i+Xl) Gi Sip''+^' + ■ 
?2asy = 52P^^+MSlP^l + ■■■ 
?3asy = Ssp^^ + uSiXiP^^'^ H 



A+e A+e+2 
a+e a+e 



and the value of the set of constants {Si,i = 1,2,3} entering the regular solution at 
large p corresponds to the value of {gi(0)} adjusted at the origin and vice versa. 

As in any eigenvalue problem, the global normalization of the eigenvector may be 
chosen at will so that, fixing gi (0) = 1 for instance, allows one to determine discrete 
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values of A. Positive values give the critical exponents, negative values are subcritical 
(or correction-to-scaling) exponents. 

The peculiar value A = 0, if present, is associated to the vestige of the 
reparametrization invariance [HZlEl- Indeed this zero eigenvalue is associated to the 
redundant operator that generates the line of equivalent fixed points in the complete 
theory jSHlEZI- Conversely, if one considers together the fixed point equations (jSHTUjl 
with the eigenvalue equations (fT^H?I|) in which A is set equal to zero (and the condition 
gi (0) = 1 is maintained), then the condition (fTTjl may be abandoned and V2 (0) adjusted 
so as to get a common solution to the set of six coupled ODE. Then, the resulting value 
of 1] nececessarily coincides with ?7°p* as defined by equation (fTH|) and the resulting value 
of V2(0) gives Zq^^. Though the number of differential equations has increased twofold 
this procedure of determining the optimized fixed point is the most efficient one when 
parameters (like and some other ones, see following section) have to be varied. 

2.4- The free parameters 

In order to perform an actual numerical study of the set of second order ODE described 
in the preceding section, I make the following choice ^: 



Following the terminology of reference jSH]; the free parameter b is redundant and 
is intended to be used to optimize the numerical results of the derivative expansion. 

The introduction of b is linked to the general property of reparametrization 
invariance which is broken by the present derivative expansion. The normalization 
ip [0) = 1 is chosen in order to distinguish the effect of simply changing the global 
normalization of the field which induces a line of equivalent fixed points (at fixed b). 
That line is customarily associated to the arbitrariness of Zq the value of Z (0) (= V2 (0)), 
i.e. the coefficient of the kinetic term in S [(f)'^] (see preceding section). Changing the 
value of b in the complete theory would induce new (equivalent) lines of equivalent fixed 
points. 

Though they are part of the same invariance, the two free parameters b and Zq 
have effects of different nature in the derivative expansion. As a global constant of 
normalization, one can expect that, at a given order of the derivative expansion, the 
variation of Zq will still reveal a vestige of the invariance of the exact theory (see 
preceding section). On the contrary, the effect of b spreads out over different orders of 
the derivative expansion. Consequently, one expects to observe a progressive restoration 
of the redundant character of b as the order of the expansion increases. Regarding 

+ In reference IIHI . I considered one additionnal parameter within the cutoff function P (named a), 
however the equations were invariant in the change b ^ b/a so that one of the two parameters was 
unnecessary. 




(22) 



(23) 
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the extremely low order considered here (O (9^)), one must not expect too much from 
varying b (see section 1^7^ . 

Notice that the cutoff function P (g^) is essentially a regulator of the integrals (jHI) 
generated by the derivative expansion. Besides the arbitrary choice in the decreasing at 
large q, the other sources of arbitrariness of the cutoff function may be included within 
the arbitrary function ip{q^). This is why the choice fj22|) does not involve any free 
parameter. 

3. Numerical study and results 

There are two different methods for numerically studying systems of coupled nonlinear 
ODE as those described above: the shooting and the relaxation methods (see, for 
example reference jlD])- Because it is the easiest to implement, only the shooting method 
is considered here (though it is less numerically stable than the relaxation method). 

3.1. The shooting method 

Considering an initial point pi where known conditions (initial conditions) are imposed 
and trial values are given to the remaining integration constants, one integrates the 
ODE system up to a second point p2 where the required conditions are checked. Using a 
Newton-Raphson algorithm, one iterates the test until the latter conditions are satisfied 
within a given accuracy. 

In the present study, the two points pi and p2 are either the origin p = and a 
large value pasy (shooting from the origin) or the reverse (shooting to the origin). 

In principle, shooting to the origin is technically better adapted to the present 
study. For example, in the case of the fixed point equations (ISHTUj) with t] fixed, one 
starts from pasy with trial values for the three constants {Gj} and the initial values of the 
three functions and of their first derivatives {vj, v^} defined by fll4flT?)|l . After integration 
up to the origin, one checks whether the three conditions (fTTHT^ are fulfilled or not. 
The system is balanced and the values of the other interesting parameters are simply 
by-products of the adjustment. For example, the value of Zq associated to the arbitrarily 
fixed value of t] is simply read as the value that V2 (0) takes on after achievement of the 
adjustment. If one wants instead to determine rj {Zq), then r] has to be considered as 
a trial parameter like {Gj} and the supplementary condition ()17p with Zq an arbitrary 
fixed number, must be fulfilled at the origin. 

Unfortunately, equations (jSHTUj) are singular at p = so that it is impossible to 
shoot to the origin. In reference |^ (where the leading order LPA was studied for small 
values of A^), the difficulty was circumvented by shooting to a point close to the origin. 
On the contrary, I have chosen to shoot from the origin because it is easy to control the 
equations starting from that point. 

Starting from the origin implies that, for a given value of ?7, the adjustable 
parameters no longer are the {Gi} but the initial values {vj (0)}, the initial values 
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Table 1. Critical exponent estimates (for d = 3) from the O {N) Wilson-Polchinski 
exact RG equation expanded up to O (9^) in the derivative expansion for two values 
of the free parameter b. 
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of the first derivatives follow from (jllflT^ . At the point pasy only three conditions 
are needed in order to balance the number of adjustable parameters. Thus one must 
eliminate the {Gi} from the six equations ()14flTm) to obtain three anonymous conditions. 
Consequently the precise knowledge of the asymptotic behavior of the regular fixed point 
solution is not necessary. 

It appears that imposing conditions like {v-' (pasy) = 0} is sufficient to determine 
the solution of interest. I first adjust the trial parameters so as to reach some not too 
large a value of pasy, then I increase pasy until the desired accuracy is reached on the 
trial parameters. 

3.2. Numerical results 

At the leading order LPA, the function vi is considered alone. Only one second order 
ODE defines the fixed point (twice this number for the eigenvalue problem). The study 
has already been done in reference [H] for d = 3 and = 1 to 4. I have extended the 
results (for o? = 3) to larger values of A^. The study is simple because: 

(i) the reparametrization invariance is automatically satisfied since, by construction, 
the coefficient of the kinetic term is supposed to be constant. 

(ii) the equation (0) with V2 = V3 = 0, does not depend on any free parameter like b. 

Consequently there is no ambiguity: to any value of A^, corresponds a unique value 
of u and u while ri=0. As already mentioned in reference j^, it is noteworthy that 
the values I obtain (see also [H]) for the critical and subcritical exponents (without 
optimization since there is no free parameters at hand) agree, "to all published digits" , 
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Table 2. Same as tabled (two other values of b). 
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Table 3. Compared to LPA on the Ihs of the arrows, the order 9^ (rhs of the arrows) 
induces a splitting into two of the subcritical exponents of degree higher than uj (values 
for d = 3, = 1 and 6 = 0.11). The LPA values may be found in gS]. 
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with the "optimized" values obtained in reference [12] from a study of an exact RG 
equation of the second family. Hence, for the numerical results I have obtained at order 
LPA for I/, uj and other subleading critical exponents uOn-, the reader is referred to table 
1 of reference 

At order (9^, once the arbitrariness of Zq has been removed via the definition of 
^opt^ the values of the critical exponents still depend on h in such a way that, for a given 
value of A^, it is impossible to define "preferred" values. For example rj (from now on, 
Tj stands for 77°^*) depends on h almost linearly. 

Due to the spreading out over several orders of the effects of 6, the study of the 
convergence of the derivative expansion relies rather on the consideration of several 
orders. Now the order is still too low to allow an appreciation of the convergence of 
the derivative expansion. Instead of presently producing one "best" estimate for each 
critical exponent at a given value of A^, it is preferable to maintain the freedom of h in 
order to better emphasize the early beginnings of some criteria of convergence if any 
(see section . 

Tables HI and El display the estimates of 77, z/ and ui as obtained for four values of b 
and A^ varying from 1 to 20 (while d = ?>). 
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For = 1, the results of reference [THj (preferred value of ?7°p* for b = 0.11) is 
obtained again. In this latter work I had proposed a criterium of choice of the value 
of b which gave a preferred value of r^. It was based on the idea of Golner of a 
global minimization of the magnitude of the function Z (0) (i.e., V2 (0)). However the 
extension of this criterion to iV > 1 is not easy to implement because the function Z {(f)) 
has been split into two parts {Z{p) and Y{p)). 

The further discussion of these results is left to section 13.31 

Table El displays a comparison between the LPA and 0(5^) results for d = N = 1 
and b = 0.11 of the series of subcritical exponents Un (with ui = u). Compared to 
LPA, the order O {d"^) has increased the potential number of subcritical exponents (due 
to the supplementary terms in the truncated action considered). For instance, if one 
adopts the dimensional analysis of perturbation theory, then at order LPA a term like 0^ 
contributing to S[(f)] is associated with the subcritical u!2 while a 0^-term generates 003. 
But at order O (9^), new terms involving two derivatives contribute to S[(f)], and a term 
proportionnal to 0^(90)^ induces a priori an in-between correction exponent. Despite 
important differences between the classical dimensions of the couples 0^ and 0^(90)^, 
table El clearly shows that, for A^ = 1, the order induces a splitting of the LPA values 
of the subcritical exponents Un for n > 1 {u is, in the perturbative approach, associated 
to the unique 0^ coupling and is not, fortunately for the perturbative framework, subject 
to this splitting). This simple splitting is presumably not preserved for other values of 
A^ due to the supplementary contribution of Y{(j)). 

3.3. Comparison with other studies, discussions and conclusion 

Figure^shows the evolution of 77 with A^ from different works. The results obtained from 
the resummation of six j29j and seven jHOI orders of the perturbation field theory serves 
the purpose of standards (other accurate estimates of the critical exponents, especially 
for A^ < 4, exist in the litterature, for a review see [Sj). One sees that the present study 
yields generally small values of rj compared to the standards except for small A^ and 
for the highest values of b. The evolution of rj with A^ is smoother than in the work 
of Morris and Turner j2B| but, as in this latter work, the non-monotonic behavior of 
the standards (responsible for the maximum of 77 about A^ = 2 or 3) is not reproduced. 
Instead, the results of Gersdorff and Wetterich j2Z| are better. The present results are 
however not so bad if one keeps in mind that the first estimate of rj in the derivative 
expansion is given at order d^. In particular, figure shows also recent estimates from 
the resummation of three orders of the perturbative series using an efficient method |44j . 
One sees that the present estimates withstand the comparison (except the monotonic 
evolution with A^). Notice also, for A^ fixed, the monotonic evolution of rj with b already 
mentioned. 

The results for the critical exponent u are more interesting to discuss because the 
order d"^ provides its second estimate. Figure El shows the evolution with A^ of at order 
compared to the results at order LPA (obtained in the present work), the standards 
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[221 iSOl are reproduced also. Again, one observes that the results for large N are not as 
good as for small values. However, for these latter values, one clearly sees that there is 
a range of values of b where the two estimates at orders LPA and flank the standards 
and another range where the two present estimates are on the same side (with respect to 
the standards). This is a phenomenon often observed in convergent series the elements 
of which depend on a free parameter (like b) but the resumed series does not: on varying 
the free parameter, one may observe monotonic or alternate approaches to the limit. 
These features may be used to determine error bars. Presently one additionnal order 
would be necessary to propose such error bars. Figure |21 shows also that, when 
increases the dependence of on 6 becomes non-monotonic. This is interesting since 
such extrema may indicate a vestige of the primarily independence on b of the exact RG 
equation. Why this effect does not occur at small values of A^ is not explained. Once 
more, at least one supplementary order of the derivative expansion would be necessary 
to understand this point. 

Figure El shows the results for u coming from i2£| and j2Zl compared to the present 
results for b = 0.11 and the standards j^niEOl- One observes that the present results 
are globally better than in [SE! and that again the estimates of Gersdorff and Wetterich 
fI7\ are excellent (for small A^ the points almost coincide with the standards). 

Figure |31 shows the present results for u at order LPA and d"^ compared to the 
standards [SHUSOl- This figure is the matching piece to figure El and the same kinds of 
remarks stand: monotonic and alternate approaches to the standards at fixed A^ exist 
as well as non-monotonic dependences on b. The magnitude of these effects are larger 
than for u and the accuracy is worse, but this is expected for a subleading eigenvalue: 
the accuracy decreases as the order of the eigenvalue increases. 

Figure El shows the results for uj coming from [20] compared to the present results 
for b = 0.11 and the standards fIE\ I30j. One observes that the present results are much 
better than in |26j . 

One may regret that Gersdorff and Wetterich |2I], who obtained excellent values 
for u and rj, had not estimated uj. As already said, this study [2ZI does not follow the 
conventional approach defined in the introduction. In particular nothing is said on the 
way the reparametrization invariance is accounted for. In fact, instead of leaving free the 
value of Z{0) = Zq to get a function i]{Zo), the procedure followed in |2I| was to attach 
the determination of r] to the minimum of the potential. This condition fixes Zq and the 
arbitrariness carried by the reparametrization invariance is implicitly removed this way. 
Because Morris and Turner have considered an equation of the same family and have 
obtained disappointing results, I think that this particular way of choosing ?7°p* could 
be the main reason for the excellent estimates of the critical exponents obtained in |27j . 
It will be interesting to adapt it to the study of the Wilson-Polchinski RG equation. 

To conclude, the derivative expansion at order 0{d^) already displays a tendency to 
converge. This must be confirmed by considering the next order which is in progress |32j . 
The study of reference which for A^ = 1 follows the procedure of reference |2Z! and 
the optimization process of reference jSHI is very encouraging. I think that the Wilson- 
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Polchinski RG equation, which is the simplest exact RG equation, is better adapted to 
the estimation of the critical exponents. Further studies should be undertaken with a 
view to better determine the status of the derivative expansion. 
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Figure 1. Critical exponent rj as function of {d ~ 3). Open and full circles 
represent the standards, they were obtained from the resummation of six (O from 
|29p or seven (• from PHj) orders of the perturbation field theory. Recent estimates 
of 1331 from the same perturbative series at order three only are also represented (A). 
The other estimates are from the exact RG equation expanded up to O (9^) in the 
derivative expansion: full triangle up from |27| . □ from |26| . The points linked by 
straight lines are from the present work: ■ for 6 = 0.03, full triangle down for 5 = 0.11, 
+ for 6 = 0.17011, * for b=0.25611. See text for discussion. 
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Figure 2. Critical exponent i' as function of N {d — 3). Open and full circles 
represent the standards as in figure ^ The other estimates are from the present work: 
□ from order LPA, the points linked by straight lines are from O (9^) in the derivative 
expansion: ■ for b — 0.03, full triangle down for b = 0.11, + for 5 = 0.17011, * for 
b=0. 25611. See text for discussion. 
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Figure 3. Critical exponent v as function oi N {d ~ 3). Open and full circles 
represent the standards as in figure ^ The other estimates are from the exact RG 
equation expanded up to O (9^) in the derivative expansion: full triangle up from 
|27j . □ from The estimates of the present study for b = 0.11 is reported (full 

triangles down linked by straight lines, see figure [5] for the other values of 6). See text 
for discussion. 
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Figure 4. Subcritical exponent oj as function of (d = 3). Open and full circles 
represent the standards as in figuren]but the opencircles (O ) are from |2S|. The other 
estimates are from the present work: □ from order LPA, the points linked by straight 
lines are from O (9^) in the derivative expansion: ■ for b = 0.03, full triangle down 
for b = 0.11, + for 6 = 0.17011, * for b=0.25611. See text for discussion. 
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Figure 5. Subcritical exponent oj as function of {d ~ 3). Open and full circles 
represent the standards as in figure 21 The other estimates are from the exact RG 
equation expanded up to O (9^) in the derivative expansion: □ from The estimates 
of the present study for b = 0.11 is reported (full triangles down linked by straight 
lines, see figure 0] for the other values of b). See text for discussion. 



